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SUMMARY 


OBJECTIVE 

The  objective  of  this  work  is  to  determine  the  performance  of  the  maximum  (or  peak) 
of  estimates  of  the  magnitude-squared  coherence  (MSC)  on  a  passive  ambiguity  surface  as  a 
detection  statistic. 

RESULTS 

The  method  chosen  to  measure  the  objective  was  the  receiver  operating  characteristic 
(ROC)  curves  which  plot  the  probability  of  detection  versus  the  false  alarm  rate.  To  obtain 
the  ROC  curves,  the  probability  distribution  function  of  the  maximum  value  of  the  MSC 
estimates  for  noise  only,  along  with  noise  and  signal  present,  are  developed.  These  probability 
distribution  and  density  functions  are  then  used  to  establish  a  Neyman-Pearson  hypotheses 
test  for  signal  detection.  The  signal  detection  results  are  presented  as  ROC  curves  which  de¬ 
pend  on  the  size  of  the  surface  over  which  the  maximum  is  taken,  the  true  coherence,  and  the 
number  of  degrees  of  freedom  of  the  individual  estimates  of  the  MSC.  These  results  are  then 
extended  to  multiple  coherence. 
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1.  INTRODUCTION 


< 


The  determination  of  the  probability  distribution  and  density  functions  of  the  maxi¬ 
mum  (or  peak)  of  estimates  (of  the  magnitude-squared  coherence)  on  a  passive  ambiguity 
surface  is  the  main  focus  of  this  report.  This  peak  value  of  the  estimate  of  the  magnitude- 
squared  coherence  (MSC)  serves  as  a  “data  reduction”  aid  in  processing  the  information  con¬ 
tained  in  the  passive  ambiguity  surface  (PAS).  We  obtain  the  probability  distribution  and 
density  functions  in  terms  of  the  size  of  the  surface  over  which  the  maximum  is  taken,  the 
true  coherence,  and  the  number  of  degrees  of  freedom  of  each  individual  MSC  estimate. 

These  probability  distribution  and  density  functions  are  then  used  to  characterize  the  signal 
detection  performance  of  the  maximum  MSC  estimate. 

In  Section  2  we  present  a  review  and  necessary  background  on  the  MSC  estimate  for 
a  single  point  on  a  PAS.  The  probability  distribution  function  of  the  maximum  value  of  the 
MSC  estimates  for  noise  only,  along  with  noise  and  signal  present,  are  developed  in  Section  3. 
These  probability  distribution  and  density  functions  are  then  used  to  establish  a  Neyman- 
Pearson  type  hypotheses  test  for  signal  detection  in  Section  4.  The  signal  detection  results 
are  presented  as  ROC  curves  which  depend  on  the  size  of  the  surface  over  which  the  maximum 
is  taken,  the  true  coherence,  and  the  number  of  degrees  of  freedom  of  the  individual  estimates. 
In  Section  5  we  extend  the  results  of  Sections  3  and  4  to  the  case  of  the  magnitude-squared 
multiple  coherence.  Section  6  concludes  with  some  final  remarks  regarding  the  relationships 
between  the  parameters  of  the  PAS  and  the  detection  performance. 


2.  BACKGROUND 

Let  X|(t)  and  X-,(t)  be  real,  zero-mean  Gaussian  random  processes  which  are  jointly 
wide-sense  stationary.  The  correlation  functions  Rjj.lt)  are  defined  by 

Rjkm  =  H{Xjtt  +  s)Xk(s)} 

for  j.  k  =  1.2.  and  the  associated  power  spectral  densities  0jj.(  f)  are  assumed  to  exist  with 


0jk(f)  = 


Rjk(  t)  exp  (  -  i  2ir  ft)  dt. 


for  j.  k  =  1.2.  The  magnitude-squared  coherence  (MSC)  function  is  defined  as  the  magnitude- 
squared  of  the  cross-spectral  density  0  j  -df)  divided  by  the  product  of  the  power  spectral 
densities  0|  |<f)0si(f).  i.e..  the  MSC  function  is  defined  by 


10  j  ->(  t )  I” 

0  |  |  (  t  )  0  -)  -»(  t  ) 


(  I) 


I 


» 


when  0jj(f)  >  0.  and  by  0  otherwise.  The  MSC  function  y~(f)  measures  the  proportion  of  the 
power  of  Xj,  at  frequency  f,  attributable  to  the  linear  regression  of  X j  on  X->  J 1 1 . 

An  estimate  of  7  (f)  is  formed  from  N  independent  (i.e.,  no  overlap)  segments  of 
data  given  by 


N 

53  F,  (f;  k)  F2*  (f;  k) 


53  IF,  (f:  k)  1“  53  IF:<f:k)|2 
k=l  k=  1 


(2) 


where  Fjlf;  k)  is  the  Fourier  coefficient  at  the  frequency  f  from  the  k1*1  discrete  Fourier 
transform  (DFT)  sample  of  the  process  X,(t).  We  note  that  F:  (f;  k)  is  a  complex  Gaussian 
random  variable. 

A  passive  ambiguity  surface  (PAS)  is  formed  by  comoensating  the  MSC  estimate  for 
a  range  of  time  delays  and  doppler  shifts.  Specifically,  to  account  for  a  time  delay  of  r  and 
a  doppler  correction  of  6  radians,  the  estimate  of  the  MSC  becomes 


53  Fj  (f:  k)  e"ik0  F2*  (f;k-r) 

k=l _ 

N  N 

23  IF,  (F.k)|2  £  |F2(f;k-r)|2 
k=  1  k=l 


(3) 


The  passive  ambiguity  surface  consists  of  a  "grid"  of  points  over  which  the  MSC  estimate  is 
computed  for  At(j  distinct  values  of  time  delay  T  =  {tj . TAtJ  anc*  *or  ^d  distinct  values 

of  doppler  ©  =  j  0  j . }  .  Thus,  the  PAS  contains  a  total  of  A  =  At(,  *  Aj  points,  each 

of  which  is  an  estimate  of  the  MSC.  Also,  the  spacing  of  the  time  delays  and  doppler  shifts 
is  assumed  to  give  independent  estimates  of  the  magnitude-squared  coherence.  Again  we 
note  that  F,  (f;  k)  e-'k^  and  F2*  (f;  k-r)  are  both  complex  Gaussian  random  variables. 

For  the  sake  of  completeness  we  review  the  distribution  and  density  functions  for 
the  estimate  of  the  MSC  y~  (r.  6 )  at  a  single  point  of  the  PAS.  The  frequency  f  is  assumed 
to  be  fixed,  and  thus,  the  dependency  of  y-  (f:  r.  0)  on  f  is  suppressed.  There  are  two  cases 
to  consider: 

I  a)  Signal  absent.  In  this  case  the  true  value  of  the  MSC  is  assumed  to  be  zero. 

The  probability  distribution  function  of  (2)  is  given  by  12) 


j 

< 


r 


F^(x  |N,72  =  0)  =  1  -(1  -x)N_1  (4) 

7 

for  0  <  x  <  1,  and  the  probability  density  function  is 

f^(x  |N,72  =  0)  =  (N-  I )  ( I  -x)N~2  (5) 

7“ 

for  0  <  x  <  1 ,  and  N  is  the  number  of  independent  data  segments.  These  follow  from  the 
fact  that  the  F;  (f;  k)  in  (2)  are  complex  Gaussian  random  variables  [2] . 

(b)  Signal  present.  In  this  case  the  true  value  of  the  MSC  is  assumed  to  be  nonzero. 
The  distribution  function  of  (2)  for  nonzero  true  coherence  is 


and  the  density  function  is 

->  i  N  v  -i  7  1-2N 

f  ^7  (x  I N,  7“)  =  (N  -  1 )  ( 1  -  7")  ( 1  -  x)N  _  ( 1  -  X7~) 

r 


,F ,  ( 1  -N.  1  -N;  l;x72) 


(7) 


where  tF|  (a.  b;  c:  Z)  is  the  Gaussian  hypergeometric  function  defined  by  [3] 


7F  |  (a.  b:  c:  Z) 


E(a)^(b)^ 

~ ^  ~ 

t=0 


and  (a )g  is  Pochhammer’s  symbol  (a)^>  -  T  (a  +  8)/T  (a),  i.e., 
(a>0  =  1 

< a =  a  (a  +  1  > ...  (a  +  t  -  1 ).  t  =  I.  2 . 


Note  that  when  either  a  or  b  is  a  negative  integer,  then  the  sFj  function  is  simply  a  poly¬ 
nomial  of  finite  degree.  We  also  note  that  the  distribution  and  density  functions  given  in  (4) 


r 


♦ 


through  (7)  remain  valid  for  the  estimate  given  in  (3)  since  the  Fourier  coefficients  are,  as 
noted  after  (3),  complex  Gaussian  random  variables. 


3.  DISTRIBUTION  OF  THE  MAXIMUM 


In  this  section  we  derive  the  distribution  function  of  the  maximum  of  the  estimate 
of  the  magnitude  squared  coherence  on  a  passive  ambiguity  surface.  The  estimates  y^  ( r ,  9) 
are  assumed  to  be  independent  of  each  other  and  the  total  size  of  the  surface  is  assumed  to 
contain  A  points.  Theorem  1  was  originally  obtained  in  [4] . 

Theorem  1.  Let  =  max  \  y-  (r,  6)  :  r  e  T.  0  e  @  |  be  the  maximum  over  A  inde¬ 
pendent  estimates  of  the  (noise  only)  MSC  where  the  true  coherence  is  assumed  to  be  zero 
<7“  =  0)  for  all  values  of  r  e  T  and  6  e  0  over  which  the  maximum  is  taken.  Then,  the 
distribution  function  of  is  given  by 


F7  (x  I  A,  N,  y~  =  0) 
A 


X  <0. 

0  <  x  <  1 . 
1  <  x , 


and  the  density  function  is 


(8) 


f  7  (X 


A.  N.  y~  =  0)  =  A  (N  -  IX  1  -x) 


n-: 


exp 


j^-Atl  -x)N_1  j,  (9) 


for  0  <  x  <  1 .  and  0  otherwise. 

Proof.  The  set  of  estimates  J  y-  (r.  0)f  are  independent  and  identically  distributed 
with  the  common  distribution  function  (4).  By  [5.  p.  36] ,  Z^  has  a  type  3  asymptotic 
distribution  given  by  (8)  since  1  -  Fix)  behaves  like  9  (xq  -  x)a  with  j3  =  1,xq=  l.and 
a  -  N  -  I .  That  the  density  function  is  given  by  ( 9)  follows  immediately  by  differentiating 

(8).  D 

We  note  that  if  the  distribution  function  (8)  for  Z^  is  expanded  into  a  Taylor 
series,  then  the  first  two  terms  of  the  series  expansion  gives  the  approximation 


F7  (x  |  A.  N.  7“  =  0)  —  l-A(l-x)^-^ 

A 

which  reduces  to  the  distribution  function  of  a  single  point  (4)  on  the  PAS  for  A  =  1. 
Figure  1  presents  graphs  of  ( 5)  and  (9)  for  various  values  of  A  when  the  true  coherence  is 
zero  ( i.e..  noise  only). 
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Figure  1. 


Density  functi  nsf7  (x 
A 


A.  N.  7-  =  0)  for  Z A  =  max  7“,  the  maxi¬ 


mum  of  tire  estimates  o  the  MSC'  with  noise  only  present.  (Surface  size  =  A; 
degrees  of  freedom  =  N;  'rue  MSC  =  7-.) 


We  now  turn  to  the  ease  when  a  signal  is  present,  that  is.  when  one  estimate  is 
based  on  data  witli  a  true  MSC  of  7-  =£  0  and  all  of  the  other  estimates  on  the  surface  have 
only  noise  present.  Thus,  one  estimate  of  the  MSC  has  a  distribution  function  given  by  (6). 
while  the  remaining  A-l  estimates  on  the  surface  have  a  distribution  function  given  by  (4). 
We  now  present 

Theorem  2.  Let  =  max  j  7-  (r,  6)  :  (r.  0)  e  T  X  ©|  be  the  maximum  over  A 

independent  estimates  of  7-  <t,  0)  where  one  estimate  is  distributed  according  to  (6)  and 
the  remaining  A-l  estimates  are  distributed  according  to  (4).  Then,  the  distribution 
function  of  is 

Fz  (x  |  A.  N.  7")  =  exp  j  -  ( A  -  1 )  (  1  -  x)N_1 1  x 
~A  L  J 


->f  |  (  -  k.  1  -  N:  1 :  X7“). 


(10) 


5 


tor  0  <  x  <  ! .  1  in-  density  function  is  given  by 


/  :  \  N 

fz  (x  1  A.  N,  7~)  =  (N-  1 )  (-  7  v  |  ( I  -  x )N -- 

~A  \1  ~  xy~  J 


exp  |  -  ( A  -  I )  ( 1  -x^]j,  1  -xy2)1^  ->F ,  (1  — N,  1  — N;  1 :  x72) 


N-: 

+  (  a  - 1 1  x  y 

k=0 


(tT^t)  :F!«-k.l-N;l:x72)|  . 


(11) 


for  0  <  x  <  1 . 

^ro°t-  Let  Tq  and  0q  be  the  location  on  the  surface  where  the  true  coherence  y~ 
is  nonzero.  Then 


Z 

~A 


max  1 7-  (r.  0)  :  (r.  0)  e  T  X  0 
max  {7"  lr0.  0q)  .  ZA.,j- 


where 


ZA_|  =  max  ^72  (r.  0)  :  (t.  0)  e  T  X  0  \  J(r0,  0q)|  ^ 

and  T  X  ©  \  1 1  rQ.  0q)  {  is  the  set  containing  A  -  1  points  which  excludes  the  point 
(tq.  0 q I •  Thus  the  distribution  function  of  ZA  is  given  by 

F7  (x  |  A.  N .  7“  I  =  Pr|z  <x} 

Ta  r  <  ~A  ' 

=  Pr{T“,r0'^0><x  ZA-l 

=  T">(x  |  N.7“)  F7  t x  |  A  -  1 .  N.  7-  =  0) 

7"  A-l 

which  follows  from  the  fact  that  7 2  (Tq.  0q)  is  independent  of  the  other  estimates 


|  7-  (r,  d )  :  (7-.  f>)  e  T  X  0  \  |(Tq.  }  •  The  distribution  function  (x  |  N.  7-)  is  given 

by  (4)  and  the  distribution  function  (x  I  A  -  1 ,  N.  7-  =  0)  is  given  by  (8)  of 

A-*  1 

Theorem  1 .  The  theorem  now  follows.  □ 

We  note  that  when  A  =  1.  the  extreme  value  distribution  function  of  given  in  ( 10) 
reduces  to  the  distribution  function  of  a  single  point  (with  a  signal  present)  on  the  PAS  as 
given  in  (6).  Plots  of  the  density  functions  (9)  for  a  signal  present  are  given  in  Figures  2  and 
3  for  various  surface  sizes  and  different  numbers  of  degrees  of  freedom.  Note  that  as  the  maxi¬ 
mum  is  taken  over  more  estimates  on  the  surface  the  density  function  "shifts”  to  the  right,  as 
would  be  expected.  Also  note  that  when  Figure  2  is  compared  to  Figure  3.  the  density  functions 
for  an  increased  number  of  degrees  of  freedom  N  (Figure  3)  have  less  “separation.”  Compar¬ 
ing  Figures  1  and  2  shows  the  change  between  noise  only  and  noise  plus  signal  present. 


Figure  2.  Density  functions  f7  (x  |  A.  N.  7“)  for  ZA  =  max  7“,  the  maximum 
~A  ~A 

of  the  estimates  of  the  MSC'  when  a  signal  is  present.  (Surface  size  =  A:  degrees 
of  freedom  =  N;  true  MSC  =  7“.  I 
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Figure  3.  Density  functions  f7  (x  I  A.  N,  y~)  for  Z.  =  max  7“,  the  maximum 
~A 


of  the  estimates  of  the  MSC  when  a  signal  is  present.  (Surface  size  =  A;  degrees 
of  freedom  =  N;  true  MSC  =  7-.) 


4.  APPLICATIONS  TO  DETECTION  STATISTICS 

In  this  section  we  apply  the  results  of  Theorems  1  and  2  to  the  problem  of  signal  de¬ 
tection  on  a  passive  ambiguity  surface  using  the  maximum  (or  peak)  value  on  the  surface  as 
a  detection  statistic.  We  formulate  a  Neyman-Pearson  hypotheses  test  and  obtain  expressions 
for  the  probability  of  false  alarm  and  probability  of  detection  in  terms  of  the  number  of 
degrees  of  freedom  N  ( i.o..  the  number  of  independent  segments  of  data),  the  size  of  the 
passive  ambiguity  surface  A.  and  the  true  coherence  y~. 

Let  the  observation  of  the  maximum  be 
=  maxjrlr.flldr.flltTXe) 

that  is,  is  the  maximum  value  of  the  estimates  of  the  MSC  on  a  surface  containing  A  points. 
The  hypothesis  Hq  of  the  Neyman-Pearson  test  is  the  case  when  no  signal  is  present,  which  is 
characterized  by  a  true  coherence  of  7-  =  0  for  each  estimate  on  a  PAS  of  size  A.  Under  Hq 


8 


hL 


the  observation  f «  =  ZA  is  described  by  the  density  function  (7)  given  in  Theorem  1.  The 
hypothesis  H  j  is  the  case  when  a  signal  is  present,  i.e.,  the  true  coherence  7^  is  nonzero  for 
one  particular  pair  of  (tq,  0q)  and  is  zero  for  the  other  A-l  points  on  the  surface.  In  this 
case  the  observ  ation  fA  =  ZA  and  has  the  density  function  (9)  of  Theorem  2.  Thus,  the 
hypotheses  test  is 


H0  :  ?A  "  ZA  ’ 

Hi  :  fA  =  5a- 


P0<x>  =  fZ  A  <x  I A'  N'  7"  =  0)  ; 

A 

p  I  ( x )  =  f^  (x  I  A.  N,  7-  #  0) . 


The  probability  of  false  alarm  and  probability  of  detection  are  now  readily  calculated  as 


%  =  pfa,xo  )=J  PQ<x)dx 


fZa(x0IA,n.7-  =  0) 


^d  "  PD,x0!  J"  pjtxldx 


=  1  “Fza(x0IA.N.7-) 

where  the  distribution  functions  F^  and  F^  are  given  in  (ft)  of  Theorem  1  and.  respectively. 

in  (8)  of  Theorem  2.  We  note  that  in  Neyman-Pearson  type  detection  the  probability  of  false 
alarm  is  set  at  some  predetermined  value.  Solving  for  the  threshold  xq  in  terms  of  the  proba¬ 
bility  of  false  alarm  Qq  gives 


<0  =  I  -  1  /A)  t!n  ( 1  -  Qq)J 


For  the  sake  of  comparison  we  also  determine  the  detection  performance  of  the  MSC  estimate 
for  a  single  point  [  ft  ]  1 7 1 .  In  this  case  the  hypotheses  are 


H0  :  ?  =  7: • 
h.  :  r  =  72- 


pn(x)  =  f^  (x  |  N.  7“  =  0)  ; 
7~ 

p  1  ( x )  =  (x  |  N.  7-  ^  0)  . 
7“ 
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where  (x  I  N.  7-  =  0)  is  given  by  ( 5)  and  i  ->  (x  |  N,  7-)  is  given  by  ( 7).  The  probability 

7“  7“ 

ot  false  alarm  and  probability  of  detection  are  easily  calculated  as  above. 

The  comparisons  between  the  detection  performance  of  the  “peak”  (or  maximum) 
value  of  the  estimates  of  the  MSC  on  a  PAS  and  a  single  estimate  of  the  MSC  are  given  in 
Figures  4  through  6.  In  these  figures  receiver  operating  characteristic  (ROC)  curves  depict 
the  detection  performance.  Figure  4  shows  that  the  detection  performance  of  the  maximum 
decreases  as  the  size  of  the  surface  increases  (with  the  degrees  of  freedom  N  and  true  coherence 
7-  both  fixed).  Figure  5  indicates  the  expected  result  that  the  probability  of  detection  in¬ 
creases  as  the  true  coherence  increases.  In  Figure  6  the  increase  in  the  probability  of  detection 
when  the  number  of  degrees  of  freedom  also  is  increased  is  indicated.  We  note  that  in  the 
above  figures  the  ROC  curves  for  A  =  1  can  be  obtained  from  [6] . 

Since  the  ROC  curves  presented  in  Figures  4  through  6  have  a  linear  scale,  the  detec¬ 
tion  performance  for  small  values  of  the  probability  of  false  alarm  is  difficult  to  determine. 
Thus,  in  Figures  7  through  9  the  probability  of  detection  is  plotted  against  10  logjQ  (7-), 
where  7-  is  the  true  coherence  and  the  probability  of  false  alarm  is  fixed.  Figure  7  again 
shows  that  the  detection  performance  of  the  maximum  degrades  as  the  size  of  the  surface 
increases.  In  Figure  8  we  see  that  increasing  the  probability  of  false  alarm  also  increases  the 
probability  of  detection.  Finally,  the  increase  in  the  probability  of  detection  as  the  number 
of.degrees  of  freedom  increases  is  indicated  in  Figure  9. 


Figure  4.  ROC  curves  tor  the  detection  statistic  =  max  y~.  (Surface  size  =  A; 
degrees  of  freedom  =  N;  true  MSC  =  7“.) 
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PROB.  DETECTION 


PR08  FALSE  ALARM 


Figure  5.  ROC  curve  for  the  detection  statistic  f  ^  =  max  y~.  (Surface  si/e  =  A. 
degrees  of  freedom  =  N:  true  MSC  =  7"  ) 


PROB  DETECTION 


A  1 


PROB  FALSE  ALARM 


f  inure  (i  K(K  curves  lor  the  detection  statistic  f  ^ 
decrees  o|  treedoni  =  N.  true  V1SC  =  •>'  r 


max  7* 


( Surface  si/e  = 


A. 


i: 


f  igure  4.  Detection  curves  lor  f  A  =  max  y~.  (Surface  si/e  =  A;  degrees  of  free¬ 
dom  =  V.  Qq  =  probability  of  false  alarm.) 


5.  I  XTENSION  TO  MULTIPLE  COHERENCE 

The  estimate  of  the  magnitude-squared  multiple  coherence  has  also  been  considered 
as  a  detection  statistic  |  7| .  In  this  section  we  present  a  brief  review  of  multiple  coherence 
and  derive  the  extreme  value  statistics  of  the  maximum  taken  over  estimates  of  the  magni¬ 
tude-squared  multiple  coherence.  Ihese  results  are  then  used  to  obtain  the  detection  per¬ 
formance  of  the  estimate  of  the  magnitude-squared  multiple  coherence. 

We  let  p  designate  the  number  of  source  channels  receiv  ing  measurements.  In 
particular,  the  case  of  p  =  2  reduces  to  the  magnitude-squared  coherence  (MSC)  studied  in 

T 

Sections  2  and  3.  let  X(t)  =  (X|tt) . X^ltH1  be  a  real,  zero  mean,  (iaussian  vector  ran¬ 

dom  process  which  is  widc-sensc  stationary.  The  correlation  matrix  is  defined  by 

Kin=[[l  V'*'' V'l]  '  [V1] 

where  fRik"!  is  the  p  by  p  matrix  with  elements  Rj|.(t).  I  he  associated  power  spectral 
densities  are  assumed  to  exist  with 
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4 


Ojk"'» 


oo 


Rjl^  (ii  exp  ( —  i  2 7r  1 1 >  tit 


tor  !  <  j.  k  <  p.  and  0 <  l  >  is  the  spectral  density  matrix 

0( n  .[ojL,nJ 

I  he  magnitude-squared  multiple  coherence  ( MSMCt  is  defined  in  tenns  of  the  power  spectral 
ciensi t x  matrix  0 1  t'l  and  its  inverse  0~  ’(  f>  where  the  elements  of  are  designated  by 

o~’*f»  =  [[Crll'(l,]J 


I  lie  MSMC  /a*!  .  |  i  |  ]  p  of  the  jtf’  channel  is  defined  by 

I . j- 1 .  i+ 1 . P"‘  =  '  -  1  (12» 

for  i  =  I .  ...  p.  We  note  that  the  MSMC  can  be  obtained  inductive!)  in  terms  of  the  pairwise 
complex  coherence  (not  magnitude-squared  coherence)  as  indicated  in  J  8 1 .  We  also  note 
that  when  p  =  2.  the  MSMC  p ^ t )  (or  y x.  |(  ft )  reduces  to  the  MSC  7- ( fl  defined  in  (  1 1  and 


M  j.  s(  f)  =  ja  s.  |  (  ft  =  7"i  1 1 


7  .  ■> 

for  notational  convenience  we  designate  the  MSMC  as  ( or  p  j-  (ft  J  and  suppress  the  de¬ 
pendence  on  j.  p.  and  f  (or.  respective!)  .  p>. 

I  he  MSMC  estimate  ^ "( 1 1  is  lormed  trom  an  estimate  ot  the  spectral  density  matrix 
CM  ft  =  and  its  inverse  ( t  ij  ~  uIPki  l  he  MSMC  estimate  is  given  bv 


Pj-i  fl  =1-1 


p)((f>  $"(l  ij 


When  a  signal  is  absent,  i.e..  the  true  MSMC  is  zero,  the  distribution  ot  the  MSMC  estimate 
‘""’V 

M  ~  is  given  bv  f  2  | 


I  7  1  v  1 p  S. g"  :  0 1 


(  \  -  I  +  p  1 


V  =  0 


1  1  - v,N-p+f+! 


1  I  vi 


■ 


s 


and  the  density  (unction  is 


f  <y(x  I  p.  N,  p2  ~  0)  = 
V- 


r  <N) _ 

r ip-  i) r <n -p  +  n 


XP':  (  1  -  x)N”P  . 


For  the  case  of  a  signal  present,  the  distribution  is  given  by  [2]  |9] 


F">  ( x  |  p.  N,  p2)  = 


(p-2  +  k)! 

(P-2)!  k! 


-  x 


X/J- 


-,F|  (-  k,  p  -  1  -  N:p  -  1 ;  xp-) 


and  the  density  function  is 


t>>(x  |  p.  N,  p2) 
M- 


r < n ) ( i  -n~)  xp  -d  -xr  p 

r  (p  -  l )  r  <n  -  p  +  i  >  <  i  -xp2). 


tF,  (p  -  1  -  N.  p  -  I  -  N;  p  -  I ;  xp“)  . 


(14) 


(15) 


(16) 


The  following  two  theorems  for  the  magnitude-squared  multiple  coherence  extend 
Theorems  1  and  2: 

Theorem  1 Let  pZ^  =  max  |j)2  (t.  6)  :  (t,  0)  e  T  X  ©j  be  the  maximum  over  A 
independent  estimates  of  the  MSMC  and  assume  that  the  true  multiple  coherence  for  each 
estimate  is  0  (i.e.,  g*-  =  0).  Then  the  distribution  function  of  pZ^  is 


FpZA  (x  |  A.  p,  N.  p-  =  0) 


x  <0. 
0<x<  1  , 
1  <x. 


(17) 


and  the  density  function  is 
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'aiZa(x  I  A.  p,N,m:  =  Q) 


^  )(  N-p+1  Ml  -  x)^  P  exp 


A  p-':  V1-^ 


(18) 


for  0  <  x  <  ! .  and  0  otherwise. 

Proof.  As  in  Theorem  1 ,  /aZ A  has  a  type  3  asymptotic  distribution  [5]  since 


1  -  1  p,  N.  /a-  =  0) 

M- 


“C"  0° ' X,N_P^' +" 


( 1  -  x)N-P+1 


as  x  -*■  1.  The  theorem  now  follows  from  the  independence  of  the  A  estimates  for  the 
MSMC.  □ 

Again  we  note  that  Theorem  1  is  a  special  case  of  Theorem  1  ’  when  p  =  2.  We 
now  proceed  to  the  case  of  a  signal  present  in  one  of  the  estimates: 

Theorem  2'.  Let  pZA  =  max  j  ju~  (r.  6)  :  (r.  d)  e  T  X  ©  [  be  the  maximum  over  A 
independent  estimates  of  the  MSMC  where  one  estimate  is  distributed  according  to  ( 1 5) 
and  the  remaining  A-l  estimates  are  distributed  according  to  (13).  Then  the  distribution 
function  of  juZA  is 


ff2a 


(x  |  A.  p,  N. m~) 


exp 


(A-i,r;‘ ),i-x)N-p+1 


/  ,  ■>  \ 

xP-1 

n-r  \ 

-XM2/ 

N 


N-p 


k=0 


(p  -  2  +  k)l 
(p-  2)!  k! 


1  -x 
I  -  xm2 


->F |  ( —  k .  p  —  I  -N:p:x/a“) 


(19) 


for  0  <  x  <  I .  The  density  function  is 
ffiZA  (x  1 A- p- N-  ^2) 


F(N)xP-2(l -x*N'P  /  1-m2 


N 


r  (  p  -  i » r  (  n  -  p  +  i )  \j_  XJU  2 


exp 


(A-nr'VxtN-p+i 
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|  ( 1  -  '~N-)F|<p-l-N.p-l-N;p-l:  xp-) 


N-p 

+ 1  a  - 1 1  x  ( i  -  x)p_“ 

k=0 


(p-  2  +  k)! 
(P  -  2)!  k! 


k 


X ->F|  (- k.  p  -  1  -  N:  p  -  1 :  Xju-)  |  (20) 

tor  0  <  x  ^  F 

Proof.  Proceeding  as  in  the  proof  of  Theorem  2.  we  obtain 

F.y  (x  |  A.  p.  N.  JU-)  =  F^(x  |p.N.g-)F  7  (x  |  A  -  1 .  p.  N,  jlT  =  0) 

p-  H  A-J 

and  the  theorem  now  follows  directly  from  ( 1 7)  of  Theorem  I '  and  (15).  D 

Plots  of  the  density  functions  f^  <x  I  A,  p.  N.  p~  =  0)  and  f^zA  (x  I  A.  P-  N.  p~) 
are  given  in  Figures  10  through  14.  Figure  10  shows  plots  of  the  density  function  ( 18)  when 
only  noise  is  present.  (See  Figure  1  for  the  case  p  =  2.)  In  Figure  1 1,  the  effect  of  increasing 
the  number  of  source  channels  p  is  shown.  The  effect  of  increasing  the  surface  size  A  is 
shown  in  Figure  12:  also,  the  effect  of  changing  the  number  of  source  channels  p  is  contrasted 
between  Figure  1  2a  and  1  2b.  The  effect  of  changing  the  true  multiple  coherence  p-  is  indi¬ 
cated  in  Figure  15.  Finally.  Figure  14  indicates  the  result  of  changing  the  number  of  degrees 
of  freedom. 

As  in  the  case  of  the  MSC  estimates  which  form  a  passive  ambiguity  surface.  Neyman- 
Pearson  type  signal  defection  hypotheses  can  be  formulated  by  using  the  maximum  value  of 
the  estimate  of  magnitude-squared  multiple  coherence  as  a  detection  statistic.  In  this  case, 
the  measurement 

p>"A  =  max  j p-  ( t.  0  :  (r,  0 )  e  T  X  0  | 
has  either  the  distribution  under  Hq  defined  by 


»0  :  ^ A  =  pZA  P0,x'  =  *nZA(X  1  A.P-N.M2  =  °) 
or.  under  11  ] .  when  a  signal  is  present 

H,  :  pfA  =  MZA.  p  |  < x )  =  1'm2a(x  1  A.p.N.p2). 
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DENSITY  FUNCTIONS 
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Figure  10.  Density  functions  f  £  (x  I  A.  p.  N.  =  0)  for  =  max  p~ .  'he 
maximum  of  the  estimates  of  the  MSMC  with  noise  only  present.  (Surface  size 
-  A:  degrees  of  freedom  =  N;  number  of  channels  =  p;  true  MSMC  =  p~.) 
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DENSITY  FUNCTIONS 


the  estimates  ot'the  MSMC  when  a  signal  is  present.  (Surface  size  =  A:  degrees  of 
freedom  =  N;  number  of  channels  =  p;  true  MSMC  =  /j“.i 


:o 


DENSITY  FUNCTIONS 


Figure  1  2a.  Density  functions  for  pZ^  (Z^).  the  maximum  of  the  estimate  of 
the  MSMC'  (MSC)  when  a  signal  is  present.  Since  p  =  2.  fi~  =  y~  and  /xZA  =  Za. 
(Surface  si/e  =  A;  degrees  of  freedom  =  N:  number  of  channels  =  p:  true  MSMC 
=  p2  I 


x 


Figure  1  2b.  Density  functions  for  jttZA.  the  maximum  of  the  estimate  of  the 
MSMC  when  a  signal  is  present.  (Surface  si/e  =  A;  degrees  of  freedom  =  N; 
number  of  channels  =  p;  true  MSMC  =  jr~.) 
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DENSITY  FUNCTIONS  £  o  DENSITY  FUNCTIONS 


Figure  13.  Density  functions  f^A  (x  I  A.  p,  N,  p")  for  pZA,  the  maximum  of 

the  MSMC  estimates  when  a  signal  is  present.  (Surface  size  =  A:  degrees  of  free¬ 
dom  =  S';  number  of  channels  =  p;  true  MSMC  =  fJ~.) 


Figure  1 4.  Density  functions  f^z^  (x  I  A.  p.  N,  p“)  for  pZA,  the  maximum  of 

the  estimates  of  the  MSMC  when  a  signal  is  present.  (Surface  size  =  A;  degrees 
of  freedom  =  N:  number  of  channels  =  p;  true  MSMC  =  p“.) 

n 


where  f^  and  {]jtZ^  are  *?iven  < 1 8)  anti  ( -0).  respectively.  The  probability  of  false 
alarm  and  probability  of  detection,  for  the  maximum  estimate  of  the  MSMC,  are  given  by 


O0  =  I  -  FfXzA(x0  1  A-  P'  N-J“2  =  0) 
^d  =  |-FgZAl!t0|A'P'N''‘2» 


where  F  's  given  by  (  I  7)  and  F  ^  is  given  by  (  19).  The  probabilities  of  false  alarm  and 
A  ' — A 

detection  for  a  single  estimate  of  the  MSMC  are  also  easily  calculated  as 


Q0  =  1  -  F^>  (x0  |  p.  N.  p-  =  0) 

0(1  =  1  -  F-^>(xq  1  p.  N.  p-) 

u  P - 

where  F^  and  F  >  are  given  bv  (  13),  and  respectively.  (15). 

P-  P- 

The  ROC  curves  for  the  peak  value  of  the  MSMC  are  presented  in  Figures  1 5 
through  1 8.  The  curves  for  A  =  1  can  be  obtained  from  [  7] .  In  Figure  1 5  the  effects  of 
increasing  the  number  of  source  channels  p  is  shown  to  decrease  the  probability  of  detec¬ 
tion.  (See  the  discussion  in  Section  6.)  Also,  the  effects  of  increasing  the  number  of  estimates 
A  over  which  the  maximum  is  taken,  i.e.,  increasing  the  surface  size,  are  again  shown  to  de¬ 
crease  the  probability  of  detection:  this  is  clearly  demonstrated  in  Figure  16.  In  Figure  15. 
the  effect  of  increasing  the  true  multiple  coherence  p-  is  shown  to  increase  the  probability 
of  detection.  The  effect  of  increasing  N  ( the  number  of  degrees  of  freedom)  is  shown  in 
Figure  16  to  increase  the  probability  of  detection. 

Figures  19  through  22  are  graphs  of  the  probability  of  detection  versus  10  log]Q 
(MSMC)  when  the  probability  of  false  alarm  is  fixed.  Figure  19  shows  that  the  effect  of  in¬ 
creasing  the  number  of  source  channels  p  results  in  a  decrease  in  the  probability  of  detec¬ 
tion  (again,  see  the  discussion  in  Section  6).  The  effect  of  increasing  the  surface  size  A  is 
shown  in  Figure  20  -  the  probability  of  detection  decreases.  Figure  21  shows  the  expected 
result  that  the  probability  of  detection  decreases  as  the  probability  of  false  alarm  decreases. 
Finally,  we  see  in  Figure  22  that  increasing  the  number  of  degrees  of  freedom  increases  the 
probability  of  detection. 

We  conclude  by  addressing  the  fact  that  the  probability  of  detection  decreases  as 
the  number  of  source  channels  p  increases  under  the  assumption  that  the  numerical  value 
of  the  multiple  coherence  remains  fixed.  This  is.  perhaps,  difficult  to  understand  since  the 
basic  definition  of  multiple  coherence  changes  as  the  number  of  source  channels  changes. 

As  indicated  in  [  7] ,  the  signal-to-noise  ratio  (SNR)  in  the  individual  channels  must  decrease 
in  order  that  the  "multiple  coherence”  stays  constant  when  the  number  of  source  channels 
increases  (under  the  assumption  of  equal  SNR  in  all  channels).  However,  this  is  only  a  partial 
explanation. 
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Figure  I  5.  ROC  curves  for  the  magnitude-squared  multiple  coherence  detection 
statistic  pf  ^  =  max  n~.  (Surface  size  =  A;  degrees  of  freedom  =  N;  number  of 
channels  =  p;  true  MSMC  =  n~.) 


To  better  understand  multiple  coherence  we  compare  the  following  two  cases  that 
are  based  on  the  assumption  that  the  SNR  in  the  individual  channels  are  equal.  The  equal 
SNR  in  each  channel  implies  that  the  magnitude-squared  multiple  coherence  (yup  for  p 
channels)  and  SNR  (denoted  R)  are  related  by  [  7 ] 


MP  (I  +  R)(I  +<p-  1)R)  '  U1) 

For  case  1.  let  p  =  2  and  =  0.07  (=  7^).  From  (21)  we  have  that  R  =  0.359  (=  SNR).  For 
case  2.  let  p  =  3  and  let  the  SNR  be  as  determined  above,  R  =  0.359.  Thus  the  magnitude- 
squared  multiple  coherence  is  equal  to  =  0. 1 1 .  We  now  compare  case  1  where  p  =  2  and 
Ms  =  0.07  with  case  2  where  p  =  3  and  M3  =  0. 1 1 .  The  comparison  is  given  in  Figure  23,  where 
the  expected  result  of  an  increase  in  the  probability  of  detection  occurs  when  the  SNR  in  the 
individual  channels  is  held  constant  and  the  number  of  channels  is  increased. 


24 


PROS.  FALSE  ALARM 


Figure  16.  ROC  curves  for  the  magnitude-squared  multiple  coherence  detection 
statistic  =  max  fa~.  (Surface  size  =  A:  degrees  of  freedom  =  N;  number  of 
channels  =  p;  true  MSMC  =  n~.) 


PROB.  DETECTION 


Figure  18.  ROC  curves  tor  the  magnitude-squared  multiple  coherence  detection 
statistic  gtf  A  =  max  (Surface  size  =  A:  degrees  of  freedom  =  N;  number  of 
channels  =  p;  true  MSMC  =  q~.) 


10  LOG  (MSMC) 


Figure  19.  Detection  curves  for  the  magnitude-squared  multiple  coherence  de¬ 
tection  statistic  jif A  =  max  n’r-.  (Surface  size  =  A;  degrees  of  freedom  =  N; 
probability  of  false  alarm  =  Qq;  number  of  channels  =  p.) 


Figure  20.  Detection  curvesjor  the  magnitude-squared  multiple  coherence  de¬ 
tection  statistic  pf  A  =  max  q-.  (Surface  si/e  =  A;  degrees  of  freedom  =  N;  prob¬ 
ability  of  false  alarm  =  Qq.  number  of  channels  =  p.) 


10  LOG  (MSMC) 


Figure  21 .  Detection  curves  for  the  magnitude-squared  multiple  coherence  de¬ 
tection  statistic  (if  A  =  max  p-.  (Surface  size  =  A;  degrees  of  freedom  =  N;  prob¬ 
ability  of  false  alarm  =  Qq;  number  of  channels  =  p.) 
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6.  CONCLUSIONS 


1  Ik-  performance  of  the  extremum  tor  peak)  on  a  passive  ambiguity  surface  as  a 
detection  statistic  has  been  derived  in  Section  4  for  the  case  of  magnitude-squared  coherence 
and.  m  Section  5.  tor  the  case  of  the  magnitude-squared  multiple  coherence.  The  following 
conclusions,  valid  for  both  cases,  can  be  drawn: 


ill  I  he  probability  of  detection  decreases  as  the  si/e  of  the  surface  increases.  An 
intuitive  explanation  of  this  behavior  is  that^as  the  si/e  of  the  surface  is  increased,  more  noise 
is  "added"  into  the  observation  f »  =  Max  y-  (or.  in  the  case  of  multiple  coherence. 

~  M  ax  (j~  i  1  his  is  a  direct  consequence  of  the  assumption  that  a  signal  is  present  in 
onlv  one  of  the  estimates  while  the  other  estimates  are  based  on  observations  which  contain 
noise  only  I  Inis  increasing  the  surface  si/e  A  implies  that  only  noisy  data  are  incorporated 
into  the  measurement  and  that  no  contribution  ol  signal  information  is  made.  As  a  conse¬ 
quence.  the  surface  si/e  should  be  chosen  to  be  as  small  as  possible  but  y  et  large  enough  to 
ensure  that  the  signal  is  contained  within  the  surla  e. 


A 

i^i  file  probability  ol  detection  increases  as  the  true  coherence  y~  increases  lor  as 
the  Hue  multiple  coherence  p-  increases!  I  his  is  an  expected  result  which  holds  true  m  the 
asc  o|  ,i  single  estimate  as  indicated  m  |  ~| . 


iti  I  he  probability  o!  deles  1 1  on  iik  leases  ,o  the  number  ol  degrees  ot  treedom  \ 
in.  leases  lot  as  the  numlvi  o|  discrete  I  ourier  transforms  aveiaged  to  obtain  an  estimate 
in.  leases!  \eam  tins  is  an  expec  ted  lesiilt  which  also  holds  true  tor  the  single  estimate 


i  4  i  I  he  piohabihly  ol  dele,  t ion  mci eases  as  the  number  ol  source  channels  is 
m.  leased  undei  the  issiinipl iqii  that  the  signal-lo-mnse-ralio  is  the  same  in  each  o|  the 
mdiv  uhi.il  .  h.umels  <  )n  the  othei  hand  the'  probability  ot  detec  lion  decreases  as  the  number 
oi  sour,  e  ban  nc  Is  o  m.  leased  milieu  the  assumption  that  the  multiple  .  oh  ere  nee  is  t  ixed 


hi 


7.  GLOSSARY 


ABBREVIATIONS 

MSC 

MSMC 

PAS 

ROC 

SNR 


SYMBOLS 

A 

2*1 

E'l  H  .  k  i 

F"A  (\  i  N.  >:  =  <J i 
7" 

F^A  I  x  !  \.  y-  I 
7  “ 

I  /  (  x  A.  N.  7  -  -  0 1 

7  A 

1  /  (  x  A .  \ .  7  "  i 

C\ 

A 

I  "six  ;  p.  N.  -  ()i 

■> 

I  "A  I  x  p .  \ .  n  ~  » 

1  /j/.  x  '  x  A  P-  N  U' 

'M/Xlx  APNP~' 

> 

7-111.7-111 

7  "  I  I  .  T.  II  I.  7  "  (  I  .  T.  0  t 


7  -  I  T .  H  1 .  7  -  I  T .11  I 
1.  .1-1.  1+1. 

*i/A 


Magnitude-squared  coherence 
Magnitude-squared  multiple  coherence 
Passive  ambiguity  surface 
Receiver  operating  characteristics 
Signal-to-noise  ratio 


Pochhaninter’s  symbol 

Sue  of  the  ambiguity  surface 

Gauss'  hypergeometric  function 

The  Fourier  coefficient  at  the  frequence  t  Irom  the  k'^1 
discrete  Fourier  transform  1  DF'I  1  of  the  process  x^t  1 

Distribution  (unction  ot  the  MSC  estimate  when  the  true 
coherence  is  zero  (signal  absent  1 

Distribution  function  of  the  MSC  estimate  when  the  true 
coherence  is  nonzero  (signal  present  1 

Distribution  function  ot  the  maximum  over  MSC  estimates 
when  the  true  coherence  is  zero  (signal  absent ) 

Distribution  tunction  ot  the  maximum  over  MSC  estimates 
when  the  true  coherence  is  nonzero  (signal  present  1 

Distribution  tunction  ot  the  MSMC  estimate  when  the  true 
multiple  coherence  is  zero  (signal  absent) 

Distribution  tunction  ot  the  MSMC  estimate  when  the  true 
multiple  coherence  is  nonzero  (signal  present  I 

Distribution  tunction  ot  the  maximum  over  MSMC  estimates 
when  the  true  coherence  is  zero  (signal  absent  I 

Distribution  tuiution  ot  the  maximum  over  MSMC  estimates 
when  the  true  coherence  is  nonzero  (signal  present  I 

Magnitude-squared  coherence  (and  its  estimate  1 

Magnitude-squared  coherence  corrected  for  time  delax  and 
doppler  land  its  estimate) 

Same  as  the  above  with  the  frequency  t  suppressed 

Multiple  coherence  ol  channel  1  compared  to  the  other 
p-  I  channels 

Maximum  over  A  estimates  ot  the  MSMC  when  no  signal 
is  preset)’ 
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Maximum  over  A  estimates  of  the  MSMC  when  a  signal  is 
present 

The  measurement  over  the  maximum  of  A  estimates  of  the 
MSMC 

The  MSMC  when  the  SNR  is  equal  in  all  channels 

Degrees  of  freedom  (or  the  number  of  independent  samples) 
in  the  MSC  or  MSMC  estimate 

Number  of  channels 

Cross-spectral  density  between  channel  j  and  channel  k 
(or  auto-spectral  density  when  j=k)  and  its  estimate 

Cross-spectral  density  matrix  with  elements  0j|df)  and  its 
estimate 

Probability  of  detection 
Probability  of  false  alarm 

Cross-correlation  function  between  channel  j  and  channel 
k  lor  auto-correlation  function  when  j  =  k) 

Cross-correlation  matrix  with  elements  Rjj(t) 

Set  of  time  delays  at  which  MSC  estimates  are  calculated 
Set  of  doppler  shifts  at  which  MSC  estimates  are  calculated 
Data  in  channel  i 
Vector  of  data  ( x  j  ( t  > . xp(t))T 

Maximum  over  A  estimates  of  the  MSC  when  no  signal 
is  present 

Maximum  over  A  estimates  of  the  MSC  when  a  signal  is 
present 

The  measurement  of  the  MSC 

The  measurement  over  the  maximum  of  A  estimates  of 
the  MSC 


8.  REFERENCES 


1 1 1  L.  H.  Koopmans.  On  the  coefficient  of  coherence  for  weakly  .slationarv  stochastic 
processes.  Ann.  Math.  Statist..  35  ( 1964).  PP-  532-549. 

(2]  N.  R.  Goodman.  Statistical  analysis  based  on  a  certain  multivariate  complex  Gaussian 
distribution  (an  introduction).  Ann.  Math.  Statist..  34  ( 1963).  pp.  152-177. 

13]  M.  Abramowitz  and  1.  A.  Stegun.  ed..  Handbook  of  Mathematical  functions  with 
Formulas.  Graphs,  and  Mathematical  Tables.  U.  S.  Government  Printing  Office. 
Washington.  DC.  1964. 

f 4  ]  W.  Marsh,  private  communication. 

[5]  B.  Epstein.  Elements  of  the  theory  of  extreme  values.  Technometrics.  2  ( 1960). 
pp.  27-41. 

(6)  G.  C.  Carter.  Receiver  operating  characteristics  for  a  linearly  thresholded  coherence 
estimation  detector.  IEEE  Trans.  Acoust..  Speech.  Signal  Processing.  ASSP-25  ( 1977). 
pp.  90-92. 


i  7 1  R.  Trueblood  and  D.  Alspaeh.  Multiple  coherence  as  a  detection  statistic.  Naval 
Ocean  Systems  Center.  Technical  Report  265.  July  1978. 

1 8 ]  L.  H.  Koopmans,  On  the  multivariate  analysis  of  weakly  stationary  stochastic  processes. 
Ann.  Math.  Statist..  35  ( 1964).  pp.  1765-1781. 

[9 1  R.  A.  Fisher.  The  general  sampling  distribution  of  the  multiple  correlation  coefficient. 
Proc.  Royal  Soc..  Series  A.  121  ( 1928).  pp.  654-673. 


33 


